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ABSTRACT 

We report on the detection of the secondary eclipse of the very-hot Jupiter OGLE-TR-56b from combined z'-band time 
series photometry obtained with the VLT and Magellan telescopes. We measure a flux decrement of 0.0363±0.0091% 
from the combined Magellan and VLT datasets, which indicates a blackbody brightness temperature of 27181jqy K, a 
very low albedo, and a small incident radiation redistribution factor, indicating a lack of strong winds in the planet's 
atmosphere. The measured secondary depth is consistent with thermal emission, but our precision is not sufficient to 
distinguish between a black-body emitting planet, or emission as predicted by models with strong optical absorbers such 
as TiO/VO. This is the first time that thermal emission from an extrasolar planet is detected at optical wavelengths 
and with ground-based telescopes. 
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1. Introduction 

Very hot Jupiters (VHJs) are giant planets orbiting very 
close to their host star, typically with periods of less 
than 2-3 days. Their proximity to the stars makes those 
planets potentially very hot, and it is therefore expected 
that they will emit detectable amounts of thermal ra- 
diation both at optical and near-infrared wave lengths 
(jLopez-Morales fc Seagerll2007t iFortnev et al.ll2008th 

Since the first detection of thermal emission from hot- 
Jupiters (|Deming et alj |2005( ). all detections have been 
made from space using the Spitzer telescope at wave- 
lengths longer than 3.6 //m, despite numerous ground- 
based attempts (e.g. iKnutson et ail 2007t iDeming et al.1 



120071: iSnellen fc Covind 120071: iSnellerJ l2005h . Anti-transit 
exoplanet atmospheric detections from the ground so far 
have been elusive, in contr ast to the rece n t grou nd-based 
primary transit result by Rcdficld et al.l {2008). This is 
unfortunate, as ground-based anti-transit detections can 
greatly facilitate hot Jupiter atmospheric studies, and ex- 
tend this field to the optical and near-infrared regions of 

th e spectrum. 

lLopez-Morales fc Seagerl ((20071) predicted the de- 
tectability of thermal emission of OGLE-TR-56b in z'-band 
(0.9 /urn). The authors found that if the planet had a bolo- 
metric albedo Ab — and a negligible atmospheric en- 
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ergy redistribution factor, the decrement in flux during sec- 
ondary eclipse could be as large as 0.05% of the total light of 
the system. Of the hot Jupiter planets known at that time, 
OGLE-TR-56b was the best suited candidate for detection. 

Here we present the results of the follow-up observations 
to test that prediction. §2 describes our observations and 
the data reduction and analysis are presented in §3. In §4 
we summarize and discuss our results. 

2. Observations 

OGLE-TR-56b is a faint target (V=16.56) located in a very 
crowded field. These factors, combined with the small ex- 
pected signal of the secondary eclipse, make a detection 
challenging. Stable photometric conditions and large aper- 
ture telescopes, which minimize stellar blends and also pro- 
vide high resolution time sampling of the target's light 
curves, are therefore necessary. 

We first monitored three secondary eclipse events of 
OGLE-TR-56b on June 25th, 30th and July 2nd, 2008 UT 
using the FORS2 instrument at the VLT-Antu telescope lo- 
cated at Paranal Observatory in Chile. FORS2 has a mosaic 
of two 2kx4k E2V CCDs with 15x15 fim pixels. We used 
the direct imaging mode of the instrument with a Gunn-z 
filter, a high resolution collimator and 2x2 binning, to ob- 
tain images with 0.125" /pixel and a field of view (FoV) of 
4.2' x4. 2' around the target. 

Only the night of July 2nd produced high enough qual- 
ity data to detect the secondary transit of OGLE-TR-56b. 
On that night we benefited from both stable and good con- 
ditions, with an average seeing of 0.72" varying by ±0.2" 
during the run. The observations were taken between the 
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Julian dates of 2454649.67536 and 2454649.88428, yielding 
354 images with exposure times between 10 and 30 seconds 
and 26 second readout times. Exposure times were adjusted 
during the night to keep the counts in the central pixel of 
the target PSF around half-well depth. The pointing was 
stable, with the target centered to within ±0.5 pixel in the 
CCD's x-direction and ±0.25 in the y-direction throughout 
the night. 

The observations on June 25th suffer from poor and 
highly variable seeing conditions, with seeing >1.3" during 
the night, and could not be used for this analysis. The night 
of July 30th had good seeing, between 0.6" and 0.22", but 
the sky conditions changed both rapidly and dramatically 
during the in-eclipse phases, while remaining stable out- 
of-eclipse. This affects the performance of our de-trending 
algorithms (see §3), which use the systematics found in the 
out-of-eclipse portions of the light curve to correct the data 
in-eclipse. The in- and out-of eclipse systematics proved to 
be very different and impossible to remove. 

We monitored another secondary eclipse event of 
OGLE-TR-56b on August 3rd, 2008 UT with the new 
MagIC-E2V instrument on the Magellan-Baade telescope 
located at Las Campanas Observatory in Chile. The 
MagIC-E2V is a two-amplifier CCD with 1024x1024, 
13x13 /an pixels, producing a FoV of 38" x38" and a res- 
olution of 0.037" per pixel in the lxl binning configura- 
tion. All the exposures were taken with a SDSS z 1 filter. 
The observations span the Julian dates of 2454682.48142 
to 2454682.69383 containing 253 images with 60 second in- 
tegration times and 5 second readout times. The sky condi- 
tions were stable and good, with an average seeing of 0.9" 
which only varied by ±0.1" during the run. 



3. Reduction and analysis 

Several stars, including the target, appear well isolated both 
in the VLT and Magellan frames. This renders our analysis 
by aperture photometry precise enough for this study. 

After correcting for bias and flat-fielding all images us- 
ing standard IRAF routines, we started the data analy- 
sis by performing DAOPHOT-type aperture photometry, 
recording the flux from all bright non-saturated stars in 
the FoV for a large and well sampled range of apertures. In 
the Magellan data, the sky background was estimated us- 
ing annuli centered around the apertures, and with a large 
and well sampled range of inner and outer ring radii. In 
the VLT data, we estimated the sky background using an 
iterative clipping method in the central third of the image, 
which contains the target and reference stars, in a manner 
similar to the calculation of sky background in DAOPHOT. 

The best comparison stars were selected separately in 
the Magellan and VLT datasets by identifying the most 
stable (i.e. minimum standard deviation) differential light 
curves between each star and OGLE-TR-56. The individual 
preliminary light curves of both OGLE-TR-56 and the se- 
lected comparisons presented clear systematic trends, with 
the VLT dataset more highly affected. These trends can 
be attributed to atmospheric effects, such as seeing, sky- 
background brightness and airmass variations, or instru- 
mental effects such as small changes of the location of the 
stars in the detector. Those effects can be modeled and re- 
duced using de-trending algorithms commonly used in tran- 
sit light curve analyses. 
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Fig. 1. Top: OGLE-TR-56 de-trended differential light 
curves. Open squares and filled circles show, respectively, 
the VLT and Magellan data. Bottom: VLT (squares) and 
Magellan (circles) light curves binned by a factor of 23 and 
14 respectively. The blue light curve corresponds to OGLE- 
TR-56, and the black and grey triangles to reference stars. 
Models with and without the best fit secondary transit are 
shown as horizontal lines. At phase 0.65 (red squares), we 
indicate the final values and uncertainties on the fit eclipse 
depth as well as the out-of-eclipse flux. 



A standard procedure in transit light curve analysis 
is to model systematics using only portions of each star's 
light curve at phases out-of-transit and then to apply the 
model to the in-transit data, to prevent removing the tran- 
sit signal itself. For data-sets with a large number of light 
curves, syste matic errors can be corrected by algorithms 
like SysRem (jTamuz et al.l l2005h . which measures trends 
common to many stars of the sample, removing each trend 
with a pass through the algorithm. 

Our field is in a very crowded region toward the Galactic 
center, limiting the number of useful unblended refer- 
ence stars and the performance of iterative algorithms. 
Therefore, we de-trend the light curves in two steps. First 
we modeled and removed a single initial linear trend com- 
mon to the target and reference stars using an im plemen- 
tation of the SysRem algorithm (jTamuz et al.ll206"5j ). which 
seeks to minimize the expression XX r ij — c i a j) 2 / ' a fji where 
Tij is the average-subtracted stellar magnitude for the ith 
star of the jth image, a is the uncertainty of m, Ci is 
an epoch-dependent parameter, and cij is a stellar depen- 
dent parameter. After this first pass, the VLT data has 
a per-point scatter of 0.13%, down from 0.27%. Any re- 
maining trends in the target are further reduced by per- 
forming a multiple linear regression fit with position in the 
detector, FHWM of the stars and airmass. This second de- 
trending pass further reduces the scatter of the VLT light 
curve to 0.116%, close to the Poisson limit of ~0.10%. For 
the Magellan data, the initial per-point scatter of the light 
curve was 0.16% and a final scatter of 0.098% was achieved 
after de-trending, close to the Poisson limit of ^0.08%. The 
optimal VLT and Magellan apertures, and the sky back- 
ground annulus for the Magellan data, were determined for 
the target and reference stars which minimized the stan- 
dard deviation of the light curve when using the apertures 
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Fig. 2. Normalized flux histograms for the portions of the 
combined VLT and Magellan data in-eclipse (dotted line) 
and out of eclipse (solid line). The bin width is 0.000363, 
coincident with the detected transit depth. 



in the above de-trending procedures. The final light curves 
are shown in Figure 1. 

3.1. Eclipse detection and error estimation 

The VLT and Magellan datasets combined contain 596 
points between phases 0.396 5 and 0.6403, based on the tran- 
sit ephemerides equation of iTorres et al.l (|2008T ). Given the 
photometric precisions per data point of 0.116% for the 
VLT data and 0.098% for the Magellan data, the 0.05% 
maximum expected depth of the thermal emission signa- 
ture of OGLE-TR-56b during secondary eclipse is at least 
a factor of two shallower than the errors of the individ- 
ual points. Therefore, the eclipse is difficult to discern by 
eye in the unbinned data (see Fig. 1) and further tests are 
necessary to confirm its presence. 

First, we used the k nown orbital period and stellar 
and planet ary radii from Torres et al] (|2008f ) to fit tran- 
sit models (jMandel & Agolll2002t ) with no limb darkening 
to the data from each night. The best fits were found using 
a Levenberg-Marquardt least-square algorithm, with the 
central phase of the transit set t o 0.5, consistent with a 
circular orbit (|Bouchv et al.l [2005h . and leaving the depth 
and the flux out of transit as free parameters. The results 
are a transit depth of 0.037±0.016% with a reduced \ 
of 0.903 for VLT and 0.036±0.011% with a reduced X 2 of 
0.926 for Magellan. The total errors in the eclipse depth 



are estimated as crjepth = a w/N + o r where a w is the 
scatte r per-point out-o f-transit and of describes the red 
noise dPont et al.ll2006f), estimated with the "praye r-bead" 
method (see iGillon et al.l l2007t ISouthworthl I2008D . to be 
1.1 xlO -4 in the VLT data and 4xl0~ 5 in the Magellan 
data. 

We finally investigated to what extent uncertainties in 
the system's parameters affect our eclipse depth measure- 
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Fig. 3. The top panel shows the derived model eclipse 
depth versus assumed central transit phase, with the plot 
approximately covering the eclipse duration. Phase 0.5 cor- 
responds to a circular orbit. The central phase is not well 
constrained within approximately 0.01, with similar tran- 
sit depths found over a phase range of 0.49 to 0.515, but 
the transit depth falls off rapidly for central phases outside 
that range. The best-fit value (black dot at phase 0.497) is 
indicated in the bottom contour plot, along with the value 
assuming a circular orbit (cross). The joint confidence con- 
tours of the best fit are also plotted in the bottom figure 
for the 68.3%, 95.5% and 99.7% levels. 



ments. Varying the impact parameter, inclin ation, planet- 
to-sta r radius ratio, and system scale from ITorres et alj 
(|2008f ) by lcr, the measured eclipse depth in the Magellan 
data changed by only ±2.4xl0 -5 or 0.22 ddepth- Thus, our 
result is largely independent of the system parameters 

Three other tests were performed to confirm the detec- 
tion. Taking the average of the 365 light curve data points 
at out-of-transit phases versus 174 points in-transit (only 
points near the bottom of the transit), we measure a depth 
of 0.0363 ± 0.0091%, which corresponds to a 3-4cr detection 
including red-noise. We further checked the detection in a 
manner s imilar to the reporte d detection of HD209458 at 
24^m bv lDeming et all (|2005[) . Figure 2 shows the first of 
those tests, i.e. histograms of the distribution of normalized 
flux for both the in-transit and out-of-transit phase inter- 
vals. The distribution of in-occultation points is shifted by 
0.0363% of the flux, as expected from the presence of the 
occultation. The results of the second test are shown in 
Figure 3. Letting the central phase and depth of the model 
transit vary, we confirm the depth and find the center of 
the occultation to be at phase 0.497jlQ'gJg, consistent with 



2 a near-circular orbit. 



4. Discussion 

A decrease in flux of 0.0363±0.0091% corresponds to a con- 
tribution from the planet to the total z'-band brightness of 
the system of 0.394±0.098 milli-magnitudes. The light de- 
tected could be stellar light reflected by the surface of the 
planet, but planetary thermal emission is the best explana- 
tion. Theoretical wo r k on Very-Hot Jup it er atmospheres 
(|Marlev et all 119991 : ISeager et all 120001 : ISudarskv et~aT1 
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Fig. 4. z'-band secondary eclipse depth of OGLE-TR- 
56b as a function of the reradiation factor / for albedos 
Ab = 0, 0.15 and 0.3 (left-side), and as a function of the 
planet's temperature (right-side), in the case of reflected 
light (dashed line), or from thermal emission assuming the 
planet emits as a black-body (solid line). The (red) filled 
and open squares indicate the expect ed model depth of th e 
transit with and without TiO/VO (|Hubenv et al.ll2003ft . 
Our measured depth is the black dot. 



2000) predict that these planets can reflect at most 30% 
of the stellar incident light (i.e. Ab = 0.3), if their atmo- 
spheres are covered by homogeneous pure silicate clouds at 
altitudes of one millibar or higher. As will be less than 0.3 
if the planets are cloud free or have patchy silicate or iron 
clouds so stellar radiation penetrates t hrough the cloud s 
and is absorbed by gaseous molecules (jHood et alj|2008h . 
Reflected light could make OGLE-TR-56b appear bright 
enough to be detected in the z'-band, but this requires Ab 
>> 0.3, and Very-Hot Jupiter atmospheres are too hot on 
the substellar side for silicate or iron clouds to form, imply- 
ing low albedos. The measured radiation most likely comes 
from thermal emission and the effective temperature of the 
planet has to be high. We estimate a brightness tempera- 
ture of T z i = 2718^07 K, based on the eclipse depth error 
bars, and assuming the planet emits thermally as a black- 
body (see Figure 4b). Ab — and / > 0.47 are also values 
consistent with the eclipse depth (See Figure 4a). This in- 
dicates that the stellar radiation is absorbed by the planet 
and not efficiently distributed throughout its atmosphere 
via winds, but i nstead is almost instantaneousl y re-radiated 
back to space (|L6pez-Morales fc Seagerll2007lh The black- 
body T z > of the planet would be 2260±19 K if the incident 
stellar radiation were completely redistributed and Ab = 0. 

Recent detections of planetary occultations at infrared 
wavelengths (~ 3-24 /im) by the Spitzer space telescope 
have revealed properties of the atmospheric structure of 
hot- Jupiter exoplanets, most notably thermal inversions 
(Hubeny et al. 2003; Burrows et al. 2007; Harrington et al. 
I2007t iKnutson et al l [20081 ; ISeager et aill2008h . thought to 
be caused by a strong optical absorber a t high altitudes 
(|Burrows et all l2007t iFortnev et all I2008D . Hot Jupiters 



with thermal inversions have bright near-infrared and in- 
frared emission lines where the planetary flux is much 
greater than a single temperature black-body model would 
predict. OGLE-TR-56b is hotter than any of the planets 
Spitzer has observed, and at the extreme hot and highly- 
irradiated end of those planets believed to have inversions. 
While our result in Fig. 4 reveals a d e pth c losest to the 
expected value from the iHubenv et al.l (|2003f) model with 
thermal inversion, more observations of the secondary tran- 
sit of OGLE-TR-56b are necessary to improve our current 
error bars and confirm this result. Also, the contrast be- 
tween models with and without thermal inversion is greater 
in the near-i nfrared, with predic ted occultation depths of 
1-3 mmags (jFortnev et al.l [2008h . making detections pos- 
sible with modern ground-based detectors. Detecting the 
secondary transit of this planet in the near-infrared will 
also further constrain the models. 

All secondary transit measurements to date have been 
made with Spitzer. A ground-based detection now extends 
exoplanet occultation science to shorter wavelengths, more 
accessible telescopes, and fainter targets. This is particu- 
larly important for the near future, as Spitzer is soon to 
run out of cryogens, severely limiting its capabilities. The 
eventual loss of Spitzer will leave a long time window before 
the launch of the James Webb Space Telescope, when oc- 
cultation measurements will rely heavily on ground-based 
measurements. At the same time, planet-finding missions 
like CoRoT and Kepler will find a wealth of new planets to 
study. Extending ground-based measurements to other very 
hot Jupiters will enable comparative exoplanetology, where 
the atmospheric structures of hot planets can be compared 
and contrasted, to continue investigation into important 
and relevant physical processes. 
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